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Abstract 

A two- level atom interacting with a single radiation mode is considered, 
without the rotating-wave approximation, in the strong coupling regime. It 
is shown that, in agreement with the recent results on Rabi oscillations in a 
Josephson junction (Y. Nakamura, Yu. A. Pashkin and J. S. Tsai, Phys. Rev. 
Lett. 87, 246601 (2001)), the Rabi frequency is indeed proportional to first 
kind integer order Bessel functions in the limit of a large number of photons 
and the dressed states are macroscopic quantum superposition states. Fo 
approach this problem analytically use is made of the dual Dyson series and 
the rotating-wave approximation. 
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A recent experimental finding on Josephson junctions Q has shown as 
Rabi oscillations happen in strong electromagnetic fields for the two-level 
model. One of the main results of this experiment was the proportionality of 
the Rabi frequency to first kind integer order Bessel functions with the order 
given by the photon number involved in the transition. 

The study of a two-level model in a cosine time-dependent perturbation in 
the strong coupling regime Q proved that in this case Rabi oscillations involve 
only odd order first kind Bessel functions and partially explains the results of 
Ref . . On the same ground we approached the problem of a two- level atom 
interacting with a single radiation mode in the strong coupling regime in Ref. 
|| . In this paper we want to extend the results of Ref. || by discussing the 
case of the experiment of Ref. Q with a large number of photons involved. 
We will show that theory and experiment indeed agree. 

The Hamiltonian we start with has the simple form as also given in Ref. 
Q (neglecting the tunneling contribution as not essential) 

H = ua^a + -^"°"3 + 9&i(pf + o) 

being u the frequency of the radiation mode, A the separation between the 
two levels of the atom, g the coupling between the radiation field and the 
atom, cti and 03 the Pauli matrices and a and the annihilation and creation 
operators for the radiation field. This apparently simple model is not exactly 
solvable unless the rotating-wave approximation is done, in this latter case the 
solution is exactly known and the model is then named the Jaynes-Cummings 
model 0]. The model (||) is able to describe the results of Ref. jjj when the 
coupling g is large with respect to the level separation. 

To reach our goal we apply the duality principle in perturbation theory as 
given in Ref. [||||. The idea is to consider, contrarily to small perturbation 
theory as given e.g. in Ref. @, as unperturbed Hamiltonian the term 

Hq = Loa^a + gcri(a) + a). 

The Schrodinger equation {% = 1) 

H U F (t) = ijU F {t) 

has the solution [f| 

U F (t) =Y / e- tE " t \[n; ax ])([n;a x }\\X)(X\ 

2 \ 2 

being E n = nu; - 2-, a x = and 

\[n;a x ]) = e^ a ~ a ^\n) 

a displaced number state pi, n an integer starting from zero the eigenvalue of 
the operator a) a and A = ±1 the eigenvalues of o\. These states represent the 



dressed state for the system but our analysis complies with the one given in 
Ref. H by a dual Dyson series otherwise, no Rabi oscillations can be obtained 
theoretically. 

At this point we are able to write down a dual Dyson series for the Hamil- 
tonian as 



U(t) = U F (t)Texp 



dt'H F (t') 



having put 



H F (t) = U F (t)-a 3 U F (t). 



It is important to note that in the dual Dyson series (|6|), as also happens in 
the small perturbation case, when there is a resonance between the two-level 
atom and the radiation field, perturbative terms appear that are unbounded 
in the limit t —* oo and the perturbation series is useless unless we are able 
to resum such terms, named "secularities" as in celestial mechanics, at all 
order. This can be done e.g. by renormalization group methods || but here 
we limit the complexity of the mathematical analysis by simply doing the 
rotating- wave approximation and ignoring any correction to it. 

Let us look at the Hamiltonian (^). It is easily realized that can be 
rewritten in the form 101 



H F — H' + Hi 



being 
A 



= f E e "^ L - [^2 J 0[«;«i])([«;«-i]ll 1 )(-l| + l[i;a-i]><[n;ai]|| - 1)(1| 
being L n the n-th Laguerre polynomial [10| and 



Hi 



A 



E 



-i(n—m)u)t 



{n\e-^ a - a]) \m)\[n-ai]){{m-a-i]\\l){-l\ + 
{n\e^ {a - ai) \m)\[n;a^i]){[m;ai}\\ - 



The Hamiltonian H' Q can be immediately diagonalized by the eigenstates 



\^n;cr) = -^[ff|[n;ai]>|l) + |[n;a_i]>|-l>] 



with eigenvalues 



A _2g2 (4 g z " 



being a = ±1. Then we can see that each level of the atom develop a band 
with an infinite subset of levels numbered by the integer number n. The eigen- 
states can be seen as macroscopic quantum superposition states (sometimes 



named in the literature as Schrodinger cat states) Qj. We will prove that the 
two-level atom shows Rabi oscillations between these states. To prove this 
result we look for a solution of the Schrodinger equation with the Hamiltonian 
Hp by taking 



that gives the equations for the amplitudes §§ 



A 



-i(E n , a -E ma ,)t e -i{m-n)u)t 



(m|e-T7( a - at V)y + He^ (a ~ at) 



At this stage we can apply the rotating-wave approximation. The resonance 
condition is given by 

E n ,a ~ E m y - (n - m)to = 

and two Rabi frequencies are obtained. For interband transitions (<r 7^ a') 
one has 



LO 



TZ = A] (n\ sinh 
while for intraband transitions one has 

W = A|(n| cosh 
By using the relation 
"2 5 



— [a — a') 



— (a — a') 

to 



\m) I 



|m)|. 



LO 



■(at -a) 



n ) 




n\ (2g 
][ — } ' 



ml \ lo 



(m\ exp 

with lI™ n \x) an associated Laguerre polynomial [10], it easy to show that 



2 V ml V lo J \ to z 



and 



2 V Tw. \ to 



e ^ 



|4--n) ^jpj I {1 + (_!)—"] 



and then, for interband transitions one can have Rabi oscillations only between 
states differing by an odd number and we write m — n = 2N + 1, while 
interband Rabi oscillations can happen only for states differing by an even 
number and we write in this case m — n = 2N. So, finally 



(n + 2N + 1)1 \to J 



2g\ 



2N+1 



-iii.-.-w.^ fig 2 ' 



\L 



(2N- 



lo* 



and 



2g 

(n + 2N)\ \uj 



2N 
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We can interpret this Rabi oscillations as involving an effective number of 
photons 2N + 1 and 2N respectively in the transitions. So, we can take n to 
be very large and N small or zero and this is in agreement with the experiment 
described in Ref. Q. This in turn means, in agreement with the experimental 
results, 



Km A 



J2N+1 



and 



TZ' « A 



J2N 



where use has been made of the Stirling approximation for the factorial n! 
e~ n n n \/2irn at large n and the equation [10 



x \ 2 



J a (2^x) = e"2 f -J L°(x) 

that holds in the limit of n going to infinity, both for 1Z and 1Z'. 

To complete this paper, we want to show how Rabi oscillations emerge 
from eqs.(14) in the limit of a large number of photons involved, when we start 
taking as initial state e.g. \0)\g) being a|0) = and (J^\g) = —\g). Indeed, 
one has 



\0)\g) =J2 a nA0)\^n,<j) 



being 



a n ,a(0) 



e 2^ \ — 



1 



2Vnl 



-l) n ] 



For interband resonance, with n 



m 



2N + 1, we get for E ny[a=l) 



E 



m,(a'- 



_i) = (2N + with even n and odd m and n > 



m. 



TZ 1Z 

a m -i{t) = a m ,-i(0) cos(— t) - ia n ,i(0) sin(-t) 

TZ TZ 

a n ,\{t) = a n ,i(0) cos(— t) - ia m _i(0) sin(— t) 



that can be put in explicit form by the coefficients (27) giving 



and finally 



a m -i(t) = -e 2u, 



W ) y/m\ 



-fg\ 2N+1 m\ .,n' 

cos(— -t)+i — a - — — — - sm(— -t) 



2 ' V (m + 2N + y 2 



n+2iV+l 



2JV+1 / ^ ^ 

— W i — rr cos — i)+ism — t) 



In the limit of a very large number m of photons we easily realize that the 
Rabi frequency is 



ft « A|J 27V+1 ^ 
and we have the oscillating amplitudes 



a m ,-i{t) = -e ( - - 7 =cos(— t) 



/g\ m 1 . . .ft , 
a m+ 2Ar + i,i(i) = e 5^ - -=«sm(— t) 
\w/ Vm! 2 



in agreement with the experimental results of Ref. (j]]. It is interesting to note 
that the probability to find the atom in one of the two levels is in any case 
proportional to a Poisson distribution. Similar expressions can be obtained for 
resonant intraband transitions (cr = a') and so we can have Rabi oscillations 
with Rabi frequencies being proportional to odd or even order first kind Bessel 
functions. The situation is quite different if instead of a second quantized 
radiation field we use a classical cosine field Q: In this case we can have 
Rabi frequency just with odd first kind Bessel functions. This is due to the 
disappearance of the band structure for a "classical" field and then, to the 
disappearance of the intraband transitions. 

The above computation gives a strong theoretical support to the exper- 
imental findings of Ref. Q. On a different ground, we can state that the 
description through dressed states as generalized in Ref. || is sound. The 
experiment realized by Nakamura et al., besides to be a first realization of a 
strongly perturbed two-level system by a radiation a field, can be seen as the 
realization of oscillations between macroscopic quantum superposition states. 
However, it is important to point out that decoherence is observed whose 
source is to be identified in view of a practical use of Josephson junction 
as gates for quantum computation. But, it is a fundamental result to have 
proven experimentally the existence of a two-level system in a strong coupling 
regime. 

In conclusion we have shown how the experimental findings in Ref. [jl| 
can be explained theoretically by the dual Dyson series and the generalized 
understanding of dressed states described in Ref. || . 

I am in debt with Yasunobu Nakamura for giving me a preprint of their 
paper before the publication. 
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